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Mechanism of confinement in low-dimensional organic conductors 
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Confinement-deconfinement transition in quarter-filled two-coupled chains comprising dimer- 
ization, repulsive interactions and interchain hopping has been demonstrated by applying the 
renormalization group method to the bosonized Hamiltonian. The confinement given by the 
irrelevant interchain hopping occurs with increasing umklapp scattering which is induced by the 
dimerization leading to effectively half-filling. It is shown that the transition originates in a 
competition between a charge gap and the renormalized interchain hopping. 

KEYWORDS: A. Organic compounds, D. Electronic structure, D. Electrical Properties 



§1. Introduction 

In low-dimensional organic conductors, repulsive in- 
teractions play an important role for electronic states 
with a gap or a pseudo gap. The anisotropy in elec- 
tric conductivity is enhanced by interactions since the 
induced pseudo gap around the Fermi surface of a single 
chain precludes electrons from hopping between chains 
There are several arguments as to whether or not 
the electrons are confined to a chain by the repulsive in- 
teraction. Away from half-filling, the confinement needs 
a large magnitude of the interaction even for the small 
limit of the interchain hopping || since the effect of the 
interchain hopping is much larger than that of the intra- 
chain interaction. However, in the case of half-filling, the 
electrons can be confined by the interaction with a mod- 
erate strength due to umklapp scattering which induces 
the charge gap |,|,|. 

Bechgaard salts of organic conductors, TMTSF and 
TMTTF, can be regarded as effectively half- filling due to 
dimerization ^,0. The optical experiments have shown 
the finite Drude weight for the TMTSF salts but not 
for the TMTTF salts although the correlation gap ex- 
ists in both salts ||. This indicates the transition from 
an insulating state with the electrons confined to chains 
to a metallic state with deconfined electrons when the 
correlation gap becomes larger than the interchain hop- 
ping In the present study, such a transition is 
elucidated by applying the renormalization group (RG) 
method to a model of quarter-filled two-coupled chains 
with dimerization. 

§2. Formulation 

We consider quarter-filled two-coupled chains given 

by 

« = + (-l)'*d] + h.c.) 
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where t, t±, and U denote energies for the intrachain 
hopping, the interchain hopping, dimerization and on- 



site repulsion, respectively, n* 



C jal C 3< jl - 



The quan- 



tity Cj a i denotes the annihilation operator of the elec- 
tron at the j-th site of the Z-th chain (/ =1, 2) with 
spin cr(=t,|). We use the Fourier transform, Cj a i — 
N~ 1/>2 'EkCkcri exp[ikja\ with the total number of sites N 
and the lattice constant a. First, the id-term is diagonal- 
ized to obtain two bands in the reduced zone, —it /2a < 
k < it /2a, and the lower band becomes effectively half- 
filled, which band is described with fermion operators, 
dkai, and is examined in the present study. Next, diago- 
nalizing the ix-term by a ka ± = (T^fccri + d ka2 )/V2, the 



kinetic term is written as Tij^ 



(C = ±) with e(fc,±) = -2[t 2 cos 2 ka + t 2 sin 2 ka] 1 ^ 2 ± 
2t±. Thus we have the following effective Hamiltonian 
PH |. The kinetic energy with the linearized dispersion 
around the Fermi surfaces, k-p± = fcp T t±/vp, is ex- 
pressed as 11% = Y,k,p,tr,c v F(P k - fc Fc) <4 pC rC afe P°< Wlth 
p(= +, — ) denoting right moving (left moving) electrons 
and v F = V2ta [1 - (t d /t) 2 } /[l + (U/t) 2 } 1 / 2 , in which 
the independence of the velocity is discarded. Coupling 
constants of interactions corresponding to forward scat- 
tering with the same and opposite directions (g^ and g 2 ), 
backward scattering (g\) and Umklapp scattering (g^) 
are given by gi± = g 2 ± — g±± — Ua, g% oc Ua(t^/t) and 
ffiii = 5211 = 54ll = where || and _L denote interactions 
for the same spin and opposite spin. 

Applying the bosonization method to electrons around 
the new Fermi points, we introduce Bose fields of phase 
variables, 9 p+ and 9 a+ (0c+ and 9s+) which ex- 

press fluctuations for the total (transverse) charge den- 
sity and spin density, respectively |12|| . The commuta- 
tion relation with conjugate phase variables is given by 
[6i,-\-(x),6 u i—(x')] = iTr6 Vll ,isgii(x — x') . In terms of these 
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phase variables, our Hamiltonian is given by 
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, (2) 



where w p((T ) = w F [l+(-)C//7rw F ] :L / 2 , if p(o .) = [H{-)J7a/7TOF 

= «s = W, K c = K s = 1, .g p((T ) = +(-)J7a and 
g 3 = Ua(2t d /t)/[l + (t d /t) 2 }. The quantity a is a cut- 
off of the order of lattice constant. In eq. (2), there are 
twelve nonlinear terms rewritten as 

-^rr 9 / dx cos \Fh\) vv cos V2ip v <p> , (3) 

where ip v ± — 9 V ± except for -0c+ = Sc+~ {8t±x/vp)/^/2 
and -0p+ = 8 P + — ^ I (2\/2)- The RG equations for K v — 
K u (l), t± = t±(l) and G„ £ ^ p > = G vp y v i(l) are given, up 
to the second order, as 
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where iv = v^/v-p and tj_ = t±(l)/vFCt~ 1 . J n is the n-th 
Bessel function. The quantity I is related to energy scale 
to or temperature T by I = \n(W/w) or ln(W/T) with 
W(= vpa^ 1 ) being of the order of band width. The ini- 
tial condition for the RG equations are given by -ft"„(0) = 
K v , G„p : „'p'(0) = g yp ,y' P ' /2ttv f and tx(0) = tj_/(wFa _1 ) 
where gc+,s+ = -gc-,s- = 0, 3c+,s- = -gc-,s+ = 
-Ua, g a +,c+ = ~9a+,c- = ~9a+,s+ = -9a+,s- =Ua, 
and g p+tC+ = 9j>+,c- = -9 P +,s+ = 9 P +,s- = 9z- We 
take a — 2a /n and discard the RG equations for the 
velocity v v . 

§3. Conflnement-deconfinement transition 

We calculate eqs. (4)- (10) numerically by choosing U, 
id, and t± as parameters. 

Figure 1 shows the Z-dependence ofJfp(Z), K a (l), Kq(1), 
K$(l) and t±(l) which exhibit four gaps. With increas- 
ing/, K p {l) decreases to zero forming a gap in the total 
charge fluctuation while Kc(l) increases infinity to in- 
duce a gap in the transverse charge fluctuation. Quan- 
tities K a {l) and Ks(l) decrease also to zero and lead to 
spin gaps for both the total and transverse spin fluctu- 
ations. The rapid increase of t_i_(J) comes from a fact 
that the term with Kq(1) in the r.h.s. of eq. (4) reduces 
to zero due to a factor Ji(8tj_). Figure 2 displays the 
corresponding Z-dependence of coupling constants. The 
main figure shows coupling constants for forward and 
backward scatterings with Gq+.s+ (curve (1)), — Gc+,s- 
(curve (2)), G C -,s+ (curve (3)), -G C -,S- (curve (4)), 
G (J+! c+ (curve (5)), -G CT+ , C - (curve (6)), G CT+! s+ (curve 
(7)) and — G ff+i s_ (curve (8)) while the inset shows those 
for the umklapp scattering with G p +,c+ (curve (9)), Gp+,o 
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^ G,y + ^cp — [2 — K v — Kq\ G u+t c p 



dl G " + ' Sp 



— G u+ ,s+Gc P ,s+ — G^ +! s-Gcp,s-, (8) 
[2-K v -Kl]G v+tSp 
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-G u+ ,c~Gc-,s P , (9) 




Fig. 1. The /-dependence of K p (l), K„(l), K c (l), K s (l) and t ± (l) 
for U/t = 5, t±/t = 0.1 and t d /t = 0.05. 
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Fig. 2. The /-dependence of coupling constants G UPiV i p i (I) for 
with U/t = 5, t±/t = 0.1 and t d /t = 0.05 where curves (1)-(12) 
are explained in the text. 



(curve (10)), — G p+! s+ (curve (11)) and G p+ ,s- (curve 
(12)). Coupling constants G p+ ,c- and — G p+! s+ increase 
rapidly and give rise to the trigger of relevance of the cou- 
pling constant Gq-,s+ as seen also in the one-dimensional 
chain. Note that the relevance of coupling constants 
Gc-,s+i — G CT+ .c- and —G a+y s+ is also obtained in the 
absence of umklapp scattering || . The relevant behav- 
iors found from the zero limit of K p , K a , 1/Kq and 
Ks exhibit the phase locking of 9 P +, 9 a +, 9c- and 9s+, 
which are given by \f29 p + = n/2, ^/29 a + — 0, \f29c- = 
and v / 2#s+ = 7r from relevant behaviors of curves (3), 
(6), (7), (10) and (11). Other coupling constants, which 
are expected to decrease [|), are still large due to the 
second order perturbation. The change of the sign of 
G ( j+ ! s+ m the renormalization process comes from the 
relevance of 9 a+ and 9 S+ . The effect of t± on cou- 
pling constants becomes large at low energies where the 
splitting of magnitudes becomes noticeable for the for- 
ward scattering (between curves (1) and (4) and between 
curves (2) and (3)), the backward scattering (curves (5), 
(6), (7) and (8)), and the umklapp scattering (curves (9), 
(10), (11) and (12)). 

In Fig. 3, the ^-dependence of t± is shown for the fixed 
t d /t = 0.05, t dfC /t(~ 0.082) and 0.1. The increase of t d 
leads to the suppression of t±(l). The case of td/t = 0.05 
shows the relevant behavior, which corresponds to de- 
confinement . The quantity t±(l) for td/t — 0.1 does not 
increase monotonically but decreases to zero after taking 
a maximum indicating confinement. The quantity tx(l) 
with a critical magnitude of td = td. c denotes the be- 
havior between the confinement and the deconfinement. 
For comparison, we show the dotted curve (the dashed) 
curve which is calculated for U = (U/t = 5, td/t — 
0.05 but (73 = as a special choice of parameter) where 
the analytical expression for the dotted curve is given 



Fig. 3. The /-dependence of t±(l) for U/t = 5 and t±/t = 0.1 
with the fixed t d /t = 0.05, t d ,cA(^ 0.082) and 0.1. The dotted 
(dashed) curve shows t±(l) for U = (U/t = 5, t d /t = 0.05 but 
</3 = 0). The inset denotes t d c as the function of t± for the fixed 
U/t = 3, 4, and 5. 



by t±(l) = t±c l . Note that the dashed curve is different 
from the case of U/t = 5 and td = 0. The dashed curve 
is evaluated to obtain t ^' , which denotes the interchain 
hopping renormalized only by the intrachain interaction, 
i.e., without the umklapp scattering. The effective in- 
terchain hopping t e ± '° is defined by t c ± '° = t exp[— / e ff,o] 
where t±(l e s t o)/t = 1 for 93 = 0. The quantity t e ± '°/t 
becomes unity for U = or t±/t = 1. The inset denotes 
the t^-dependence of i<j,c for the fixed U/t = 3, 4, and 5. 




Fig. 4. The /-dependence of K p (l) for U/t = 5 with the fixed 
t d /t = 0.05 and 0.1 where tj_ = for the dashed curve and 
t± =0.1 for the solid curve. The arrow denotes /a denned by 
K p (Ia) = Kp/2. The inset exhibits the charge gap: Aj D for 
t±/t = (dashed curve) and A p for tj_/t = 0.1 (solid curve). 
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Fig. 5. The t j_ -dependence of tjj with U/t = 3 (dashed curve) 
and 5 (solid curve) where t^/t = (curves (1) and (4)), 0.05 
(curves (2) and (5)) and 0.1 (curves (3) and (6)). 



The boundary between confinement (t<j > id, c ) and de- 
confinement (id < td.c) depends appreciably on U, which 
gives rise to the enhancement of the confined region on 
the plane of t± and t<j. c - The limiting form for small ij_ 
is given by t^ tC /t oc (t±/t) F with the [/-dependent F. 

The Z-dependence of K p (l) is shown by solid curve for 
the fixed id/i = 0.05 and 0.1 in Fig. 4. The charge gap 
is defined by A p = v-pa~ l exp[— Za] with K p (Ia) = K p /2 
and vpa- 1 = t(7r/V5)[l - (*d/*) 2 ]/[l + (*d/*) 2 ] 1/2 - Th c 
dashed curve, which denotes K p (l) for ij_ = 0, leads 
to the charge gap, A* D , for one-dimensional (ID) case 
in the presence of dimerization, id- The charge gap is 
suppressed slightly by the interchain hopping since K p (l) 
for the solid curve decreases slowly compared with that 
for the dashed curve, i.e., Za is increased by t±. Such a 
behavior is understood from eq. (5) with v = p, in which 
the first term of the r.h.s. becomes small due to the Bessel 
function suppressed by the large t±. The inset exhibits 
the charge gap as a function of id with the fixed U/t = 
3, 4 and 5 where the solid (dashed) curve corresponds to 
the case of t±/t = 0.1 (0). Note that the dashed curve 
is given by A p D ~ W(gz/W)W- 2K ^ @. 

In Fig. 5, the independence of t e f with U/t = 5 (U/t 
= 3) is shown by the solid (dotted) curve for id/i = 
(1), 0.05 (2) and 0.1 (3) (t d /i = (4), 0.05 (5) and 
0.1 (6)), where the effective interchain hopping t±, in- 
cluding the effect of the dimerization, is estimated by 
tf = texp[-Zeff] with i±(Z off )/i = 1 @. The curves 
(1) and (4) for small ij_/i well reproduce the analytical 
result given by if /t± oc (tj_/i) a o/(i-ao) with the In- 
dependent quantity a = (K p + K^ 1 + K a + K^ 1 - 4)/4 
JD]. Comparing the slope of curve (1) with that of curve 
(4) , one finds that the renormalization of ij_ increases by 
the intrachain interaction U. With decreasing t±/t, the 
ratio tf /t± decreases rapidly and becomes zero for t± 




tj t 

Fig. 6. The phase diagram of confinement (I) and deoonfincmcnt 
(II) on the plane of tj_ and A p for U/t = 3 (dashed curve), 
4 (dotted curve) and 5 (solid curve). The inset denotes the 
corresponding phase diagram on the plane of t e ^' and A p D . 



less than a critical value of ij_ indicated by the arrow. 

From the inset of Fig. 3 and that of Fig. 4, we ob- 
tain the phase diagram of confinement (I) and decon- 
finement (II) on the plane of t± and A p . In Fig. 6, the 
boundary between these two states is shown for U/t — 
3 (dashed curve), 4 (dotted curve) and 5 (solid curve). 
The boundary is rather straight and the [/-dependence 
becomes small compared with that of the inset of Fig. 3. 
It has been claimed previously that such a behavior in- 
dicates the competition between the charge gap and the 
interchain hopping t± However such a state- 

ment is not clear enough since the boundaries depend 
on the choice of U. This problem can be remedied by 
the following treatment. We take A p D (t c f'°) in stead 
of A p (t_L ) as the vertical (horizontal) axis on the phase 
diagram where A p D is obtained from the dashed curve of 
the inset of Fig. 4. The resultant boundaries are shown 
in the inset of Fig. 6, where the good coincidence is ob- 
tained among these three curves. Thus it is concluded 
that the boundary between confinement and deconfine- 
ment is determined by the competition between the one- 
dimensional charge gap (i.e., in the absence of the inter- 
chain hopping) and the interchain hopping renormalized 
only by the intrachain interaction (i.e., without the umk- 
lapp scattering). 

§4. Discussion 

We have examined the mechanism of confinement in 
terms of quarter-filled two-coupled chains with dimeriza- 
tion as a model of low dimensional systems. The confine- 
ment occurs when the charge gap induced by the umk- 
lapp scattering becomes larger than the interchain hop- 
ping renormalized by the intrachain hopping. It has been 
found that the ratio for U/t = 5 (U/t = 3) and 0.05 < 
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Fig. 7. The phase diagram of confinement (I) and dcconfincmcnt 
(II) on the plane of t± and A p . The main figure denotes the 
boundary for U/t = 5 and t d /t = 0.05 with the fixed nearest- 
neighbor interaction: V/t = 0, 0.5 and 1. In the inset, the dotted 
(dashed) curve denotes the boundary for the fixed t<j/t = 0.05 
(0.1) while the solid curve denotes the boundary for the another 
case of the Hubbard model at half-filling with t^/t = 0. 
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t x /t < 0.2 is given by Aj D /if '° ~ 1.8 (A* D /if -° ~ 1.9) 
and 1.1 < A p /tx ^ 1.3 (1.3 < A p /tx ^ 1.5). 

Here we discuss the effect of nearest neighbor interac- 
tion (V), for the small V, in which the commensurabil- 
ity energy may be negligibly small |l4|,[l^]. In this case, 

V does not contribute to both the backward scattering 
and the umklapp scattering even in the presence of the 
dimerization (t^). The effect of V appears in the forward 
scattering where the coupling constant is expressed as 
32" = 9a\\ = 2 Va and = 94^ = (U + 2V)a. In 
Fig. 7, the boundary of confinement-deconhncmcnt tran- 
sition with U/t = 5 is shown for V/t = 0, 0.5 and 1 where 

V enhances the confined region. 

Finally we comment on the case of two-coupled chain 
with the conventional half-filled Hubbard model where 
the phase diagram is shown in the inset of Fig. 7. Com- 
pared with those of Fig. 6, the ratio A p /ij_ ( ~ 1.7 for 
tx/t = 0.1) is slightly large and the boundary (solid 
curve) is rather straight. Such a result is compared with 
the boundaries for the fixed ta/t = 0.05 and 0.1, which 
are shown by the dotted curve and dashed curve, respec- 
tively. With increasing the dimerization, the slope of 
the boundary becomes steep and moves toward the solid 
curve since the large dimerization may lead the system 
to half-filling. 
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